Super-aging correlation function and ergodicity breaking in logarithmic potentials 
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We consider an overdamped Brownian particle moving in a confining asymptotically logarithmic 
potential, which supports a normalized Boltzmann equilibrium density. We derive analytical expres- 
sions for the two-time correlation function and the fluctuations of the time-averaged position of the 
particle for large but finite times. We characterize the occurrence of aging and nonergodic behavior 
as a function of the depth of the potential, and support our predictions with extensive Langevin 
simulations. While the Boltzmann measure is used to obtain stationary correlation functions, we 
show how the non-normalizable infinite covariant density is related to the super-aging behavior. 
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INTRODUCTION 



Ergodicity is a central concept in the theory of stochas- 
tic processes. A random variable A is said to be er- 
godic when its time average A(t) — (1/i) J*dt'A(t') 
over a single realization is equal to the ensemble average 
(A) = J dA AP(A) over many realizations of the process 
in the limit of infinitely long times Here P(A) is the 
equilibrium probability density for the random variable 
A. For ergodic variables, the width of the probability 
distribution of the random time average vanishes in the 
long-time limit, and the distribution reduces asymptoti- 
cally to a delta function centered on the ensemble aver- 
age, Q (A,t) — > S (A — (A)y For nonergodic variables, 
the time average remains a stochastic quantity even in 
the infinite time limit. A general criterion for the er- 
godicity of a process is given by the Khinchin theorem 
[2J, which asserts that a stationary process is ergodic if 
its autocorrelation function (A(t)A(to)) — > (A) 2 when 
\t — to\ — > oo. However, for processes for which a station- 
ary autocorrelation function does not exist, we cannot 
use the Khinchin theorem to predict the ergodic proper- 
ties of the process. For specific systems, there are gener- 
alizations of the Khinchin theorem to the nonstationary 
case [1, H[ , but for a general system, the ergodic proper- 
ties are not straightforward to predict. Meanwhile, for 
finite times, the time average will be a random variable 
for any process, ergodic or not. Since in all experiments 
the measurement time might be large, but is always fi- 
nite, it is essential to determine the properties of the 
distribution of the time average, as the latter cannot be 
evaluated using the equilibrium measure P{A). The dis- 
tribution of the time average has been investigated for 
continuous time random walk models Q , but no general 
theory exists. In Ref. [6[, we have provided a general 
expression for the variance of the time-averaged position 
x(t) for an overdamped Brownian particle in a binding 
field and showed that, in the special case of the logarith- 
mic potential, the ergodic hypothesis breaks down. 

In the present paper, we study in more detail the vari- 



sition of an overdamped Brownian particle moving in 
a confining potential that is asymptotically logarithmic. 
This system defines an important class of processes which 
has found widespread applications in the description of 
the dynamics of particles near a long, charged polymer 
]3 , momentum diffusion in dissipative optical lattices to- 
la ] , probe particles in one-dimensional driven fluids [16| , 
self-gravitating Brownian particles |17| , long-range inter- 
acting systems fl8l - 420j and diffusion of fractals [2j. a s 
well as the dynamics of bubbles in DNA molecul es [22 - 
25[ , of vortices [26| , and of trapped nanoparticles [27j. A 
striking characteristic of these systems is that the equilib- 
rium probability distribution possesses a power-law tail 
which may lead to diverging moments and ergodicity 
breaking [28| . Their anomalous behavior is controlled 
by a single parameter that is related to the depth of the 
potential. In the following, we derive explicit long-time 
expressions for the two-time correlation function and for 
the time-averaged variance of the position of the particle, 
and investigate in detail the diffusive and ergodic prop- 
erties of the system for large, but finite times. We show 
how a super-aging correlation function describes this sys- 
tem, even when the stationary Boltzmann distribution is 
normalizable. 

Figure [T] shows a typical trajectory of a Brownian par- 
ticle in a logarithmic potential in the nonergodic phase, 
obeying the overdamped Langevin equation 
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ance a (t) 



■(t)) — (x(t)) of the time-averaged po- 



with U(x) = (Uq/2) ln(l + x 2 ) and the fluctuating Gaus- 
sian noise F(t), (F(t)F(t')) = (2fc s r/(m 7 ))<5(t - £')• 
Here, m is the mass of the particle, 7 the friction co- 
efficient, ks the Boltzmann constant and T the tempera- 
ture. Obviously the time average x(t) = (1/t) J* dt'x(t') 
does not converge to the ensemble average (x) = even 
for long times. The reason for this behavior is the long 
excursions of the particle into the tails of the potential, 
where the slope and thus the restoring force tend to zero 
for x ^> 1. On average, these excursions get ever longer 
as time increases and thus dominate the time average 
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for all times. Our goal is to characterize the nonergodic 
behavior of this process. 

The outline of the paper is as follows. In Section [Til 
we solve the Fokkcr-Planck equation corresponding to 
Eq. ([1} for the Brownian particle by transforming it into 
a Schrodinger-like equation and using an eigenfunction 
expansion. In Section IIIII we explicitly construct the 
conditional probability density from the eigenfunction ex- 
pansion and determine its asymptotic long-time behav- 
ior, as well as those of the first and second moments of 
the position variable. We then use the conditional prob- 
ability density in Section IIVI to compute the two-time 
position correlation function, which exhibits either sta- 
tionary behavior or nonstationary, super-aging behavior, 
depending on the depth of the potential. Finally, in Sec- 
tion |V1 we evaluate the variance of the time-averaged 
position and show the existence of a threshold for the 
potential depth below which ergodicity is broken. Parts 
of our results were obtained in Ref. [6( with the help of 
a scaling ansatz. We here provide a complementary and 
more detailed derivation. 



II. FOKKER-PLANCK EQUATION FOR THE 
LOGARITHMIC POTENTIAL 



We consider an overdamped Brownian particle mov- 
ing in a symmetric (U(x) = U(—x)) potential U(x) 
that is asymptotically logarithmic for large \x\ 3> a, 
U(x) ~ Uq ln(|a;|/a); during the calculation, we will set 
a = 1. For o ^ 1, the variable x in our results should 
be replaced by x/a. An example for such a potential 
would be U(x) — [7 /21n(l + x 2 ), which is what we use 
for our numerical Langevin simulations. The dynamics 
of the particle is governed by the Fokker-Planck equation 
for the probability density W(x,t) corresponding to the 
Langevin equation (Q}, 



d_ 

di 



W(x,t) 



D d 

kWT dx 



(U'(x)W(x,t)) 



D^W(x,t) 



(2) 

where D = ksT / (rwy) is the diffusion coefficient. In 
the following, we set the mass of the particle to m = 1. 
The stationary equilibrium solution to Eq. @ is given 
by W eq (x) = exp(—U(x)/(kBT))/Z with the normaliz- 
ing partition function Z — j dxex.p(— U(x)/(kBT)). As 
we will see, the partition function Z and thus the station- 
ary probability density do not always exist for the loga- 
rithmic potential. Our goal is to evaluate the long-time 
behavior of the variance of the time-averaged position, 
x(t) = f Q dt' x(t')/t, which is given by, 



a 2 (t) = (x 2 (t)) 
1 



w> 2 

ft 

dt" I dt'C(t",t') 



dt'(x(t')} 



(3) 
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FIG. 1: (color online) Typical trajectories (gray) and their 
time average (red) for the Brownian particle moving in an 
asymptotically logarithmic potential U(x) — (7o/21n(l + x 2 ) 
(top panel, inset). The parameter a measures the ratio of the 
depth of the potential Uo and the temperature (see Eq. IjlOj l). 
here we used ksT = 0.5, 7 = 1. The top panel shows the 
nonergodic phase (Uo = 2fcsT), where the time average x(t) 
does not converge to the ensemble average (x) (which is zero, 
dashed line). The bottom panel shows the ergodic phase 
(Uo — 4fcsT); here the time average does tend to the en- 
semble average for long times. Also note the different scales 
for the a;- axes. 



For an ergodic system we expect x(t) — > (x) cq = and 
thus a 2 (t) — > 0. The two-time position correlation func- 
tion C{t, to) = (x(t)x(to)) in Eq. ([3]) can be expressed in 
the form I29II. 



C(t,t ) 



dx / dxa xxqP(x, t\xo, tn)W(xo, to) 



(4) 



where P(x, t\xo, to) is the conditional probability density. 
This expression is correct for a prior initial condition 
of a narrow diffusing packet (e.g. a Gaussian distribu- 
tion centered on the origin) and sufficiently long times 
to- For power-law initial conditions, a different behav- 
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ior is expected, see Refs. [3CJ, |31| . For most potentials 
U(x) (e.g. harmonic), W(x,t ) in Eq. (jl]) can be re- 
placed by the stationary probability density W oq (x). For 
the logarithmic potential, however, it turns out that the 
integrals in Eq. Q do not always converge and we need 
to use the explicitly time-dependent probability density 
W(x, to). In order to solve the Fokker-Planck equation 
@, we transform it into a Schrodinger equation and em- 
ploy an eigenfunction expansion. Writing the probabil- 
ity density as W(x,t) — xi x )^{ x ,t) with the function 
X(x) — exp (—U(x)/(2knT)), we obtain the following 
equation for ip(x, t) 29], 



19,, . d 2 ,, . 



Ddt 



1 



4fc|T 2 



dx 2 
U l2 {x) + 



1 



2k B T 



U"(x))l>{x,t). (5) 



Equation ([5]) has the form of an imaginary-time 
Schrodinger equation with the effective potential, 



U e g(x) 



1 



Ak\T 2 



U l2 (x) + 



f 



-JJ'\x) 



2k B T 

Its general solution is given by the expansion, 



ip(x,t) 



A 



-D\t 



where if)\ (x) are the eigenfunctions of Eq. ([5]) , 



dx 2 



U e g{x) 



4>\{x) = \->p\{x) . 



(6) 



(7) 



(8) 



Since the Hamiltonian is symmetric in x, we look for 
solutions of the eigenvalue equation ((5]) that have either 
even or odd parity [32]. To simplify the calculations we 
use the potential 



U{x) = UoH\x\)Q{\x\-l) , 



(9) 



which has the desired asymptotic form U(x) ~ Uo In 
|x| > 1, and is zero for |x| < 1. The asymptotic long-time 
behavior of the system is up to a constant factor indepen- 
dent of the potential near the origin and our calculation 
can be generalized to arbitrary, asymptotically logarith- 
mic potentials [33]. For |a;| > 1, the eigenfunctions are 
given by @, 



(10) 

hkJ-a(\kx\)) , 



ipkA x ) = A k\J\ x \ (aik J a {\kx\) + a 2k J- a (\kx\)) 

M x ) = A(a 1 \x\^- a + a 2 \x\ a +^ , 
ipk,o( x ) = B k sgn{kx)^J\x\{bi k J a {\kx\) 
with k = ±vA and 

Uo , 1 



2k R T 



(11) 



The subscript e (o) refers to the even (odd) solutions, 
J a (x) denotes the Bessel function of the first kind, A, 



A), and B k are normalization constants and Oj, aik and 
hk € {1,2}) will be determined by the boundary 
conditions at x — ±1. The spectrum consists of a sin- 
gle discrete solution for k = 0, iJjq(x), which is an even 
function of x and for a > 1 corresponds to the ground 
state, and a continuum of excited states for k > 0. The 
parameter a > 1/2, which measures the ratio of the po- 
tential depth Uq and the temperature T, will turn out to 
be the key quantity controlling the long-time behavior of 
the particle. Using the eigenstates (fTU|) . the solution of 
the Fokker-Planck equation ([2]) is 

W(x,t) = a x( x )*Po( x ) 

poo 

+ X(x) / dka k ^ e {x)e~ Dk2t 



+ X( X ) / dkb k ipkA x ) e 



-Dk 2 t 



(12) 



The coefficients a k and b k are determined by the initial 
condition. For \x\ < 1, the potential vanishes and thus 
the corresponding eigenfunctions are 

^k,e( x ) = &k cos(fca;) , 

4>o(x) = A , 

4>k,o( x ) = B k sin(kx) , (13) 

with the boundary conditions at x = 1, 

V>fc,e(l) = ^fc,e(l) ) 

U 



2k R T 



^,e(l) 



(14) 



The last term on the right-hand side stems from the dis- 
continuity of U'{x) at x = 1 (similar boundary conditions 
are obtained for the ground state and the odd eigenfunc- 
tions). From this, the coefficients for the ground state 
are easily seen to be a\ = 1 and a 2 = 0. 

The two integrals in the solution (|12j) can be regarded 
as Laplace transforms. The large t behavior of the sys- 
tem is then determined by the small k expansion of the 
eigenfunctions, according to the final value theorem 34|. 
Using the Taylor expansion of the Bessel function 35], 
(Eq. (9.1.7)), 



J a {k) 



(15) 



T{a + 1) 

and those of the sine and cosine, we find from Eqs. (|10p . 
(fT"3l and (fl4"l) the coefficients of the excited states to lead- 
ing order in fc, 

2a - 1 (}£ " 



ai k 



-r(«)- 



i 



a 2k ^ r(i 



'->2k 



T(a) 



r(l-a) 



2a - 1 



a+l 



(16) 
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The probability density W(x, t) should be properly 
normalized at all times. For the discrete ground state, 
the normalization integral reads, 



2 J dx%{x)+2j dx%(x) = l. (17) 

As it turns out, the ground state is only normalizablc for 
a > 1. The corresponding normalization constant is, 



A 



a-1 1 



2a - 1 Z 



(18) 



Thus the discrete solution is tpa{x) = (1/V^) |x|~ Q+1 / 2 , 
which is precisely the square root of the normalized sta- 
tionary solution of the Fokker-Plack equation @ for 
|x| > 1. For a < 1, there is no normalizable ground 
state and thus no stationary probability density. Since 
the excited states form a continuum, they need to be 
normalized to a delta function, 



2 / dxil) kye / {x)\l) k ^ e/o {x) 



+ 2 J dxt/i k)e/o (x)'ip k > ie / (x)=5(k-k') . (19) 

Realizing that the main contribution to the normaliza- 
tion integral ([T9]) comes from large values of x, we use 
the asymptotic expansion of the Bessel function [35| . 
(Eq. (9.2.1)), 



\ —7— cos (ka 

V nkx \ 



to obtain the expressions, 



1 



ira 
~2~ 



(20) 



4 a lk + a 2k + 2a lfc ft 2fc COS (2(j) a - §) 
k 1 



4 b lk + b lk + ^> lk b 2k cos (24> a - f ) ' 



(21) 



where <p a = an/2 + 7r/4. Since the coefficients and 
bj k are only known up to leading order in k, we further 
need to expand the normalization constants A k and B k 
to leading order in k in order to be consistent. Since 
a > 1/2, bi k wn l always dominate 62/c and we have 



Bl 



kj_ 



(22) 



However, the relative magnitude of a\ k and ai k depends 
on the value of a. For a > 1, a\ k dominates a,2k, and 
vice versa for a < 1. This leads to 



A-1 



k 1 

k 1 
4 37 



for 
for 



a > 1 
a < 1 



(23) 



The final piece of information needed to fully determine 
the solution ([12")) are the expansion coefficients and b k . 
They are related to the initial condition W(x, 0) by, 



0-k 



a = 



dxip k ,e{x) 

dxipoix) 



W(x, 0) 
W(x,0) 



dx ipk,o{x) 



W(x, 0) 



(24) 



III. CONDITIONAL PROBABILITY DENSITY 

The evaluation of the correlation function Q requires 
the computation of the conditional probability density 
P(x, t\xo, to), the solution of the Fokker-Planck equation 
([2]) with the initial condition P(x,to) — S(x — xq) [29|. 
Accordingly, it is of the form (TT2"j) with the time variable 
t replaced by the time difference 

t = t - t . (25) 

The expansion coefficients are then, from Eq. ([24]) . 

ipk,e( x o) 



00 



bk = 



x(xo) 

4>kA x o) 

x(xo) 



A 



(26) 



As a result, we have for x > 1 and xq > 1 (the behavior 
for negative x and xq follows from symmetry), 

P(x,t\x ,to) = Ax(x)ip (x) 
X(x) 



X(xo) jo 
X{x) 



X&o) Jo 



dk t/jk,e(x )ip k A x ) e Dh T 
dk ■4>k,o(xo)4 , k,o( x ) e ~ Dk2 ' r 



A 2 x 1 ~ 2a 



2 X Xri 



dk 



A\ (a lk J a (kx ) + a 2 kJ- a {kxo)) 



-Dk?T 



x (aifc J a {kx) + a 2 kJ-a(kx)) e 
+ Bl (bi k J a (kx a ) + b 2 kJ-a(kx )) 
x (b lk J a (kx) + b 2k J- a {kx)) e~ Dk2r 



(27) 



For completeness, the discussion of the cases where ei- 
ther x or Xo is smaller than 1 is given in Appendix [21 
The very first term on the right-hand side of Eq. ([2~T1) is 
the contribution of the stationary state and only appears 
for a > 1 (since A = for a < 1). The expansion of 
the product in the integral yields four terms for the even 
contributions (those containing the Oife-coefficients) and 
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four terms for the odd contributions (those containing 
the ^-coefficients). Since the leading order of the nor- 
malization constant Ak is different depending on whether 
a is larger or smaller than 1 (see Eq. (|23p). we have to 
examine twelve integrals in total (some of which are for- 
tunately the same). However, all these integrals share 
the same basic structure: 

f°° 2 

x / duu X 'i°^-» J t _ t (2uy)J„(2uz)e~ u , 
Jo 

(28) 

where fi and v are equal to ±a and both the 
M-independent factor c e / ofJil/ (z,y) and the exponent 
A e / p „ are different depending on the combination of 
/i and v. Here we have introduced the variables u = 
fc(flr) 1 / 2 , z = x/^Dt) 1 ' 2 and y = x /(4L>r) 1 / 2 . We 
can then write Eq. (j2"7) as, 



P(x, t\xo, to) — A x 



2 -2a+l 



+ E E (^) 



2-A, 



^e/o,iJ,A z 'y) ■ ( 29 ) 



It is important to note that the main contribution to the 
integrals 3 e / ^ „(z, j/) comes from small values of u as 
the exponential factor causes the integrand to vanish ex- 
ponentially for large values of u and all other functions 
increase at most as a power of u. Table U gives a sum- 
mary of all the values of /z, v and A e / M „, as well as the 
prefactors c e / 0i(tt)I/ (z, y) of the respective integrals. 



A. Probability density function 

For very large times and fixed ieo, the conditional 
probability density loses its dependence on the ini- 
tial condition and reduces to the probability density, 
P(x,t\x ,t ) ~ W(x,t) (see Eq. (gU below). In this 
limit, the variable y — xq/^Dt) 1 / 2 is small, and the 
Bessel function containing y in Eq. (|2T))) can be expanded 
using Eq. (fT5j) to get 

•^e/o,fj,,v{ z t y) — _|_ ^ c e/o,/j,^ i z i y) y^ 

POO 

x / duu K '"-^ + ^ 1 J v {2uz)e- u2 . (30) 
Jo 

Note that the variable z = a;/(4Dr) 1 / 2 is not necessar- 
ily small, since we are interested in the large x regime. 
Equation (f3"0"|) allows us to find the leading order con- 
tribution to W(x, t) at very large times without having 
to compute explicitly all the integrals. Comparing the 
values for /i and A e / Q M _„ for the different integrals, one 
finds that the leading contribution stems from the even 
terms with /i = —a and v = a for a > 1, and from the 
even terms containing fi = —a and v = —a for a < I. 
Evaluating the respective integrals using Eq. (6. 631. 1) of 



Ref. [3(| , we obtain the asymptotic behavior of the prob- 
ability density for x > 1, 



W(x,t) 



ZT{a 



1 x l-2a 



for a > 1 



r(i-a) (4£>t)°~ V- 2a e~ro^ for a < 1. 

(31) 



Here r(a,a;) is the incomplete Gamma function. Equa- 
tion (22) for a > I is the infinite covariant density (ICD) 
[3ll l33l |37| . This result depends on the specific form of 
the potential U (x) only through the partition function 
Z, as long as it is regular at the origin and has the same 
asymptotically logarithmic behavior as Eq. ©. The ICD 
is non-normalizable, but allows computing the asymp- 
totic long-time behavior of the moments of order 
q > 2a — 2 [37]. The lower order moments are finite for 
the equilibrium density W eq (x) — e^ u ^^ kBT ^ / Z and 
thus can be obtained from the latter. For a < I, the 
asymptotic form (|3ip of the probability density is nor- 
malizable and can be used to calculate all moments. The 
asymptotic behavior of the second moment is given by 



(x 2 {t)) =2 dxx 2 W(x,T) 
Jo 

a—l 



3a-6 



for a > 2 



Z r( tt )( 2 - tt ) (^) 2 " tt forl< a <2 
(I - a)4L>r for a < 1 . 



(32) 



For a > 2, the second moment tends to a constant which 
depends explicitly on the specific shape of the potential 
for x of O(l) and can be obtained from the equilibrium 
distribution. For a < 2, the second moment increases 
with time and diverges in the infinite time limit. In the 
range 1 < a < 2, we observe subdiffusive behavior, which 
depends on the shape of the potential for x of 0(1) only 
through the partition function Z. When a < 1, the dif- 
fusion is normal and asymptotically independent of the 
potential at small x. As it turns out, the ICD Eq. (l3"Tj) 
not only governs the time-dependence of the higher or- 
der moments, but is also crucial for the calculation of the 
correlation function (see section IIV[) . 



B. Conditional probability density 

In order to calculate the conditional probability den- 
sity, we consider times long enough that we can use the 
small-A: expansion of the coefficients, while not necessar- 
ily so long that we can ignore the initial condition. In 
this regime neither z = x/ (ADt) 1 ^ 2 nor y = xq/ '{ADt) 1 / 2 
appearing in Eq. (f2"5)) are small. In the case of both z 
and y (and thus x and xq) being of the same order, not 
all the integrals can be explicitly evaluated. However, we 
can estimate the term which is of leading order for large 
times by evaluating the integrals for z = y. Fortunately, 
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prefactor c e/o ^ >v {z,y) 


*^e/ o,(j,,i/ 




V 








• 1 


1 


a 


a 


even contribution, a > 1 


y a z- a+1 X < 


2 2 - 2cv r(2-a) 
r(a)(2a-l) 

2 2_2 °r(2-a) 
r(a)(2c«-l) 


2a — 1 
2a- 1 


—a 


a 








/ 2 2~2 Qr(2 _ a) \ 2 
I \ r(o)(2a-l) i 


4a -3 


—a 


— a 






✓ 


/2 2 °- 2 r(a)(2 Q -l)y 
V r(2-a) y 


—4a + 5 


a 


a 








2 2cv ~ 2 r(a)(2a-l) 


—2a + 3 
-2a + 3 


a 

— a 


—a 
a 


tven CUI11I ID ULlOIl, UL V* ± 


j/ Q 2 - a+1 X 




r(2-a) 
2 2a ~ 2 r(a)(2a-l) 








r(2-a) 








1 


i 
i 


— a 


— a 






1 




i 
i 


a 


a 






2~ 2a r(l-o,)(2a-l) 








odd contribution 


y a Z - a+1 x < 


r(a+l) 
2~ 2a r(l-c«)(2a-l) 


2a + 1 


a 


— a 










2a + 1 


— a 


a 






( 


' 2 - 2a r(l-a)(2a-l)'\ 2 


4a + 1 


—a 


—a 








v r(a+i) ) 



TABLE I: Prefactors and parameters for the integrals appearing in the conditional probability density, Eq. (|29p . 



the integrals for the leading order term can be computed 
explicitly. For both z and y at least of order unity, the 
leading order term for long times is then (see Ref. [36j . 
Eq. (6.633.2) for the computation of the integrals), 

P e {x,t\x ,t ) 



±^1-20 



1 rrl — a a 



(ADt 

exp (- x2+x » \ I fsa. 



for a > 1 



(4Dr) 



-1 — a ry,a 



x exp 



/ x 2 +zg \ J ( xxp \ 
\ ADt J 1 ~ a \2DtJ 



for a < 1 



(33) 



for the even contribution, I a (x) is the modified Bessel 
function of the first kind, and 

P o (x,t\x ,t ) 

,(4I, T )-V-.* ; exp(-^)4(|^) ,34) 

for the odd contribution. The latter is the same as the 
probability density for the Bessel process with an ab- 
sorbing boundary condition at the origin (2(| [38| . The 
Bessel process describes diffusion in a purely logarithmic 
potential (and thus a l/x-force) which is not regular at 
the origin. While the Bessel process correctly describes 
the asymptotic behavior of the odd moments, its Green's 
function is well denned only for a < 1, when the singu- 
larity on the origin is integrable. Thus for the calculation 
of the even moments (including normalization) for a > 1 , 
the regularization of the potential on the origin is vital, 
since only then do we have a finite partition function Z. 



The odd part of the probability density does not depend 
on Z at all and is identical to the Bessel process. 

If one of the variables z and y is much bigger than the 
other, the Bessel function containing this variable will 
oscillate rapidly (see e.g. Eq. (|2T))) ') for all but very small 
values of the variable of integration u. We may then ex- 
pand the Bessel function containing the other variable 
for small arguments. For the leading order even contri- 
bution, we obtain the ICD Eq. (13"TT) if z ^> y (x ^> xq), 
and 

P e (x,t\x ,t ) 

r(1 1 _ a) (4£>t)"- V- 2a e~roV for a < 1 

(35) 

if y > z (x > x). Note that Eqs. and J3S]) arc not 
limiting forms of Eq. (|33p . since the leading order term 
is different depending on the size of z and y. The lead- 
ing order odd contribution is always given by Eq. (|34l) . 
Using the odd part of the conditional probability density 
we can directly calculate the average E(xo,t) of the 
position x (for xq > 1), 

/>oo 

E(x ,r)=2 dxxP o (x,t\x Q ,t ) 



2T(a 



^ ( 4Dr)»-#.-A lFl (|, a + 1 ,^.). 

(36) 



Here iFi(a, b, x) is the confluent hypergeometric func- 
tion. For long times, 4Dt ^ Xq, the hypergeometric 



7 



function is approximately unity, and we thus find, 



E(x ,t) 



a „ 2 a 



2T(a + l) 



(37) 



For a = 1/2 (Uq = 0), we recover E(xo,t) — xo, which 
corresponds to free diffusion. For a > 1/2, on the other 
hand, we observe an algebraic dependence on the initial 
position xq. 



IV. CORRELATION FUNCTION 

Having derived the asymptotic behavior of the condi- 
tional probability density and of the probability density 
function, we can next evaluate the correlation function 
dl]). Since the probability density W(xo,t ), as given by 
Eq. (I31[) . is an even function of xo, only the odd part 
of the conditional probability density contributes to the 
integral of the correlation function: 



C(t,to) 



dx 



dx xx P {x, t\xQ,to)W{xo, t ) 

Bcsscl process ICD 

(38) 



This result is remarkably intuitive: For a > 1, the in- 
finite covariant density (ICD) gives the probability for 
finding the particle at xq, while the Bessel process de- 
scribes the relaxation of the particle's average position 
from xq toward the origin. Since we use the asymptotic 
expressions for the conditional probability density and 
the probability density, the resulting expression for the 
correlation function is also asymptotic in the sense that 
we require both r = t — to and to to be large. Inserting 
the results Eqs. ([31]) and ([34]) into Eq. ([38]i. we obtain 
for a > 1, 



C(Mo) 



-(4Dr)~ 



dx 



" ZT(a 

( x 2 + x\ \ ( xxo \ _ 

xexp {—wT M2^) r 



dx x 2 ~ a x\- a 



4Dtr 



(39) 



With the help of the variables z — xj '{ADt) 1 / 2 and y 
xq/(4:Dt) 1 / 2 introduced above, Eq. ([39]) simplifies to 



C(Mo) 



ZT(a) 



(4Dr) 



2— a 



dz 



2-ct„,2-a 



dyz 2 - a y 



/4Dt V4D 
~2 , 2 



x e 



-(z 2 +y 2 ) 



I a (2zy)T{a,y 2 - 



(40) 



For long times, t» 1, we may take the lower boundary of 
the integrals to 0, since the integrand vanishes as z 2 and 



y 2 (note that I a (x) ~ (a;/2) Q /r(a + l) when x -> 0). The 
integral over z can then be evaluated using Eq. (6.631.1) 
of Ref. [36J, together with I a (x) = (—i) a J a (ix), to give, 



C(t,t 



-(4Dr) 2 



ZT(a + l)r(a) 
J™dyy 2 e-y 2 ^ Q, a + 1, y 2 ^j T (a, y 2 



(41) 



At the same time, for a < 1, we have the integral, 
4 



C(t,*o) 



(4D T ) 2 ~ a (4Dt ) 



a-l 



r(i-a) 

ilr. I dyz 2 - a y 2 ~ a 



/4Dr V4Dt 
,2 , .,2n 



x e- (z2+a2) / a (2z 2 /)e" y * 



For large r, the latter expression reduces to 



(42) 



c ^) g r(^R (4Dr) U 

x J" dyy 2 e-y\¥ 1 Q,q + e""^ . (43) 



Using r = t — to, the structure of the results ([4"TT) and 
(1431) can be summarized in the compact form: 



' zr(a ^ )r(a) (4g(t-*o)) 2 -%(^) 
for a > 1 



C(Mo)- ^ 



r(a + lff(l-a) 4£>(t-f )ga(^) 

for a < 1 



(44) 



where we have introduced the two functions, 

Us) = J™ dy y 2 e~y 2 1 F 1 Q, a + 1, y 2 ^j T (a, y 2 s) , 
g a {s)=s 1 - a £ dyy 2 e~y 2 1 ¥ 1 Q, a + l,/) e 



(45) 



Note that the correlation function depends on the specific 
shape of the potential for |x| < 1 only via the partition 
function Z (for a > 1) or not at all (for a < 1). These 
results are thus valid for arbitrary (regular) potentials 
with the same asymptotic behavior as Eq. (O. The be- 
haviors of f a (s) and g a (s) for small and large arguments 
(corresponding to to 3> t — to and t — to 3> to) are given 
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in Appendix [Bl For t — to ^> to, we obtain, 



C(t,t ) 



v^r( a +|) 
3zr(o+i)r(a) 



2 " (4^o) 2- " for a > 1 



4r(a+l)r(l-a) 

x (^) 2 ° (4D* ) for a < 1 



(46) 

In both cases, the correlation function is nonstationary, 
its value increases with the initial time irj and decays 
as (t — to) 1 / 2 ~ a - We note that this is the same time- 
dependence as the first moment E(xq,t), Eq. ([57]). For 
to 3> t — to, on the other hand, we find, 



C(Mo) 



zr(»j(2-«) (^DtoY 
(I - a)4Dt 



for 1< a < 2 (47) 



for a < 1 



For a > 2, the correlation function is stationary in this 
limit. This is illustrated in Fig. [5] For times t— to that are 
short compared to to, we observe the stationary behavior 
Eq. (|4"T1) . For longer times, there is a transition to the 
nonstationary behavior Eq. (|46)l. though this is difficult 
to observe in our Langevin simulations. If we start out 
in the stationary state (which corresponds to to = oo), 
we have the stationary behavior Eq. (|47p at all times for 
a > 2. 

For a < 2, the asymptotic behavior of the correlation 
function is nonstationary and is dominated by the in- 
crease of the second moment (x 2 (t )), Eq. (l32|) . which 
means that using the stationary distribution W cq (x) to 
calculate the correlation function gives an infinite value. 
In this case, the correlation function can be expressed in 
a convenient way by normalizing it to the value of the 



second moment 



C(t,t Q ) 



forl<a<2 



v5 



r(a+l)r(2-a) 



7, / t-tg 



for a < 1 



with 



Us) 



„2-a 



U(s) , 



J a (s) = sg a (s) 



(48) 



(49) 



In this regime, the system exhibits aging. However, 
contrary to usual aging behavior, which is of the form 
C(t,t ) = (x 2 ) c J((t-to)/t ) I, SlIU, Eq. m shows 
that here C(t, t ) = (x 2 (t )} f ({t — t a )/t )- Since the pref- 
actor increases with time, we call this behavior super- 
aging. A similar behavior, albeit with a logarithmic time- 
dependence has been observed for Sinai's model for dif- 
fusion in a random environment [4l| . The stationary 



Z> °< 01 

1E-3 




= 2.5 

— simulation, t = 1 000 

— simulation, t = » 

— correlation function, t = 1000 
stationary correlation function 

— non-stationary asymptote 
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t-t 



FIG. 2: (color online) The correlation function C(t,to) for 
U(x) = (U /2) ln(l + x 2 ) (thus Z = 0Fr(a - 1)/T(a - 1/2)) 
and a = 2.5 (U = 1, k B T = 0.25, 7 = 1). The solid black line 
is the analytical expression (|44p . The red and green lines are 
the results of Langevin simulations for two different times, 
to = 1000 and to = 00 (i.e. starting from the equilibrium 
distribution). For to = 1000, we observe a transition from the 
stationary behavior Eq. (|47|) to the aging form Eq. (|46|l . For 
to = 00, we observe the stationary correlation function (|47[) 
at all times t. 



si 0,1 
o* 



■ 1.5 

— simulation, t Q = 300 

— simulation, t = 3000 

— normalized correlation function 



FIG. 3: (color online) The normalized correlation function 
C(t,t )/{x 2 (to)} for U(x) = (C/ /2)ln(l + x 2 ) and a = 1.5 
(Uo = 1, k B T = 0.5, 7 = 1). The solid black line is the 
analytical result (|48p . The numerical simulations for to — 300 
(red) and to = 3000 (green) perfectly confirm the super-aging 
behavior. 



correlation function for a > 2 (|47[) agrees with the re- 
sult derived using the equilibrium solution of the Fokker- 
Planck equation 9]. However, neither the super-aging 
behavior nor the long time limit Eq. (|4"6"|) can be ob- 
tained from the equilibrium distribution: both require 
the infinite covariant density (|31D . Figure [3] shows the 
normalized correlation function (|48p for different values 
of to. It clearly illustrates the super-aging behavior of 
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the correlation function. 



V. VARIANCE OF THE TIME AVERAGE 

We are now in a position to evaluate the long-time be- 
havior of the variance a 2 (t) of the time-averaged position 



a 2 (t) = (x 2 (t)) - (x(t)) 2 
We begin with the second term, 
1 



t - 1* 



t \ 2 

dt'E(x ,t') 



(50) 



(51) 



involving the integral over E(xo,t) — (x(xo,t)). Taking 
the time average of the first moment (|37p , we obtain, 



(x(f)> 



3 - 2a 4D(i - **) 



(4Di' 



(52) 



Since we only know the asymptotic long-time behavior 
of quantities like the first moment E(xo,t) and the cor- 
relation function C'(t,to) — the short time behavior will 
in general depend on the shape of the potential U (x) for 
small x — we calculate the time average starting at some 
time t* , which we assume to be large enough that the 
asymptotic expressions are valid. For large times, t^$>t*~, 
we then have, 



<x(i)> 2 



\l-2a 



for a < 



jt^xi^Dtr^Dt*)*-^ for a > § 



(53) 



The analysis of the first term, 



(x 2 (t)) 



dt" / dt'C{t".t') 



is carried out in detail in Appendix [C] We obtain 
(4D*)- 1 for a > 3 



(54) 



a 2 (t) 



with 



Crv = < 



\t)) ~ c a { (4Dt) 2 - a for 1< a < 3 (55) 
ADt for a < 1 , 



VE f °° A„ s* " f ( \ 

Zr(o+l)r(a)(4-o) JO u (s+l) 4 -° Ja\*J 

for 1 < a < 3 

3T(a)f(l-a) fo ^ S (s+1) 3 9a(s) 

for a < 1 , 



(56) 



which dominates the contribution from Eq. (|53p for long 
times and thus determines the behavior of & 2 (t). The 



prefactor c a for a > 3 depends on the specific choice 
of the potential (see below). Note that the asymptotic 
result for a 2 (t) is independent of t* as it should. 

The explicit expression (l55|) for the time-averaged vari- 
ance a 2 (t) gives important information about the ergod- 
icity of the diffusion process in a logarithmic potential: 
the process is (mean)-ergodic if and only if a 2 (t) van- 
ishes at long times An analysis of Eq. (j5"5")) reveals 
that a 2 (t) ~ t x when t — > oo, with X = 1 for a < 1, 
X = 2 — a for 1 < a < 3, and X = —1 for a > 3 (see 
also Figs. [Hand [5]). Accordingly, the process is ergodic 
for a > 2 (ksT/Uo < 1/3), while ergodicity is broken for 
a < 2 (1cbT/Uo > 1/3). The ergodicity of the process 
for a > 2 is in agreement with the Khinchin theorem 
, since in this case there exists a stationary correlation 
function Eq. (|4T)l which vanishes as time tends to infinity. 
For 2 < a < 3 (1/5 < k B T/U < 1/3), even though the 
system is ergodic, the slow decay of <r 2 (t) means that 
ergodicity is reached anomalously slowly. The noner- 
godic behavior for a < 1 is not surprising, since in this 
regime, there exists no stationary equilibrium state. For 
1 < a < 2, on the other hand, we observe broken ergodic- 
ity even though the Boltzmann equilibrium distribution 
is normalizable. This is due to the fact that the sec- 
ond moment increases with time [37| and thus diverges 
for the equilibrium distribution. The asymptotic power- 
law exponent of a 2 (t) is shown in Fig. HJ The analytic 
prediction (|55|) perfectly matches the numerical results, 
except near the points a = 1 and a = 3, where the con- 
vergence of the Langevin simulations is very slow and we 
expect transient logarithmic corrections. The asymptotic 
algebraic behavior of the time-averaged variance is fur- 
ther confirmed in the double logarithmic plot presented 
in Fig. [5] We again observe perfect agreement between 
analytics and numerics for different values of the param- 
eter a. 

For a < 3, the results in Eq. (|5"5"|) can be fully obtained 
from our asymptotic analysis and depend on the explicit 
form of the potential U(x) only through the partition 
function Z (for 1 < a < 3) or not at all (for a < 1). 
For a > 3, on the other hand, the contribution we get 
from the asymptotic behavior of the potential is of the 
same order as the one stemming from the behavior of 
the potential for x ~ 0(1). The asymptotic analysis, 
while correctly predicting the & 2 (t) ex t~ l behavior, thus 
fails in reproducing the prefactor c a . In terms of the 
asymptotic analysis, c a appears to depend on the time 
scale A (see Appendix [CJ , which determines how long 
it takes for the asymptotic results to accurately describe 
the behavior of the system and thus can only be obtained 
by comparing the asymptotic results to the exact solu- 
tion. However, c a can be obtained from the Boltzmann 
equilibrium distribution W cq (x) 6], which by definition 
depends on the shape of the potential U (x) in the whole 
space. Figure [S] shows the prefactor of a 2 (t) for the spe- 
cific choice of the potential U(x) — i7o/21n(l + x 2 ). This 
can be interpreted as a generalized diffusion coefficient 
for the time average: a 2 (t) ~ c a (ADt) x . 
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J i I i L 



0,5 1,0 1,5 2,0 2,5 3,0 3,5 

a=U /(2 k B T)+1/2 

FIG. 4: (color online) Power-law exponent of the long-time 
variance a 2 (t) ~ t , Eq. (|55|l (red solid line), and numerical 
data obtained from Langevin simulations by fitting the long- 
time behavior up to t = 4000 (black dots), as a function 
of the parameter a = Uo/(2JcbT) + 1/2, for the potential 
U(x) = ((7 /2)ln(l + x 2 ). 
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1 2 3 4 5 

a 

FIG. 6: Prefactor c a of a 2 (t) ~ c a {ADt) x , Eq. l|55]l, as a 
function of a for U(x) = ({7 /2)ln(l + a; 2 ). Note that it is 
zero at the transition from normal to subdiffusion (a — 1) 
and diverges at a = 3. Inset: Detail for a < 2. The result for 
a > 3 was taken from [jj. 



10000 



I b 




100 1000 10000 

t 

FIG. 5: (color online) Asymptotic long-time behavior of the 
variance a 2 (t) simulated for various values of the parameter 
a = U /(2k B T) + 1/2, for the potential U{x) = (U /2) ln(l + 
x 2 ). The dashed lines are the analytical predictions given by 
Eq. (|55p . The simulation data is for ksT — 1, 7 = 1. 



VI. CONCLUSION 

We have performed a detailed analysis of Brownian 
motion in an asymptotically logarithmic potential. We 
have obtained explicit expressions for the position corre- 
lation function (1441) and the variance of the time-averaged 
position (|55|) . The asymptotic time-dependence of the 
these quantities is determined by the single parameter 
a = U~o/(2kBT) + 1/2, which measures the ratio of the 
potential depth Uq and the temperature. Both diffusion 
and ergodic properties are controlled by a, and the sys- 
tem exhibits a surprising variety of different behaviors. 



The ergodic and aging properties of the system are closely 
related to the occurrence of the non-normalizable infinite 
covariant density (ICD). 

For shallow potentials, a < 1, diffusion is normal which 
implies (x 2 (t)) cx t, a nonstationary aging position cor- 
relation and broken ergodicity a 2 (t) cx t. In this case, 
there is no stationary distribution, so the time-dependent 
solution determines all moments. For medium strength 
potentials, 1 < a < 2, the system exhibits subdiffusion 
[lf|, (x 2 (t)) oc t 2 ~ a , while the correlation function is still 
aging and ergodicity is broken, d 2 {t) cx t 2 ~ a . In this 
regime, the non-normalizable ICD determines the second 
moment of the position. The ICD is also essential for 
determining the correlation function, as the Boltzmann 
equilibrium density yields an infinite value for the latter. 
We find a super-aging correlation function C(t. to), which 
behaves as C(t,t ) ~ (x 2 (t ))f((t — t )/to) with a prefac- 
tor (x 2 (to)} which increases with time. For deep poten- 
tials, 2 < a < 3, diffusion freezes, (x 2 (t)) cx constant, and 
the correlation function has a stationary limit C(t,to) ~ 
C(t — to), which can now be determined from the Boltz- 
mann equilibrium distribution. In this regime, the pro- 
cess is ergodic, although the decay of a 2 (t) cx t 2 ~ a is 
slow. Only for very deep potentials, a > 3 do we obtain 
cf 2 (t) cx t~ x , which is the behavior expected from the 
usual (e.g. harmonic) confining potentials @- 
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Appendix A: Discussion of the center part of the 
probability density and correlation function 



In the main body of the work, we have only discussed 
the tail part, that is the case x > 1 and xq > 1, of the 
probability density P(x,t\xo,to), arguing that the con- 
tribution of the other cases (x < 1 or i < 1, or both) 
to the correlation function is negligibly small. In order 
to determine the probability density for those cases, we 
have to use expression (|13l) instead of (fTO]) for the wave 
functions including the respective variables in the expan- 
sion (|27p. So, for example, for the case x > 1 and xq < 1, 
we have for Eq. (l27l) . 



„1— a poo 



P(x,t\x ,t ) ~ A'x L - Za + 2-^ 



dk 



ipo(xo) Jo 

x A\i>k,e{xQ) {ai k J a {kx) + a 2k J- a (kx)) e~ Dk2r 

+ Bli> Ko {x ) (b lk J a {kx) + b 2k J- a (kx)) e- Dk ^ 

(Al) 

We now employ the same analysis as we have used for 
the tail part (see section Hill and in particular Eq. (12T))) ') 
in order to determine the leading order for long times. 
We then use a small-A: expansion of the even and odd 
wave functions in the center, Eq. (JTT 



(A2) 



i>k,e( X o) — A k , 4>k,o{ x o) — B k kx 



Using Eq. (|A2[) (with the coefficients given by Eqs. (IT6"1) . 
0, (USD and 36], Eq. (6.631.1)), we obtain for x > 1, 
xq < 1 to leading order in r, once again dividing into 
even and odd contribution, 

P e (x,t\x ,t ) 

WR^- 2Q r(a,^) fora>l 

r(i-a) (iDT) a - 1 x 1 - 2a e-^ for a < 1 , 

2 



P o (x,t\x ,t ) 



{ADT)- l - a xx e-^ . (A3) 



The other cases can be obtained in the same way and we 
find for x < 1, xq > 1, 

P e (x,t\x ,t ) 

zTR r («> Trk) for a > 1 

- < 

^ rp L 5y (4Dr)°'- 1 e-i* for a < 1 , 
2 



P (a;,t|xo,to) 



r(a 



■(4Dr)- 1 - a M§ a e-'* 



(A4) 



and for x < 1, xo < 1, 
P e (x,t|xo,to) 



v for a > 1 

rTT ^ y (4Cr)«- 1 fora<l, 



P (a;,t|a;o,*o) 



4a 



(ADry^xxo ■ (A5) 



The even part of Eq. (|A5[) also yields the long time limit 
of the probability density P(x,t) for x < 1. Using the 
above expressions for the probability density, we can cal- 
culate the three remaining contributions to the correla- 
tion function Eq. ((31)]) for long times and a > 1, 

/OO />1 
dx / dxoxx P (x,i|xo,io)lU(£o,io) 

2^ 



3ZT(a) 



(4D(t-t Q ))*- a 



rl p OO 

4 / / dx £XoP (x,t|xo,to)W(xo,*o) 



3ZT 2 (a) 



(4D(t-«o)) 3_a fta 



t - 1 



to 



r-1 <-l 

4 / <ix / c?x a;xoP (x,t|xo,to)W / (xo,io) 



[i JO 

8a 
9ZT(a) 



(4£>(t-to))- 



(A6) 



with 



M«)= / rfyy 2 e~ a r(a,y 2 s) , (A7) 
Jo 

which is approximately constant for small s and behaves 
as s -3 / 2 for large s. For a < 1, we then have 

/OO /> 1 

dx / dxo xxoP (x 7 t\xo,to)W(xo,to) 

20F 



■(4£(t-io))2- Q (4L>t r- 1 



3r(a)r(l - a 



"1 /"OO 

4 / dx dx Q xxoP (x,t\xo,t )W(x ,t Q ) 



20? 



■(4i9(t-< ))t- Q (4Dto) Q - 1 



3r(a)r(l - a 



f - k 

~to~ 



r-1 /-l 

4 / dx / dx xx P (a;,i|xo,to)W / (a;o,io) 



8a 



9r(a)r(i - a 



■(4i9(<-t ))- i - Q (4^o) a - i , (A8) 
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X > l,xo > 1 


X > 1,X() < I 


a; < l,xo > 1 


x- < 1 , Xq < 1 


T > to 


1 


3 

t 2 


3 

r 2 


3 a 

r-h 2 


to > T 


1 


s 

t 2 5 a > 2 
r 2 t Q , a < 2 


r~i, a>2 
r 2 i Q , a < 2 


r" 3 , a>2 



TABLE II: The contributions from the different parts of the correlation function relative to the contribution from x > 1, xq > 1 



with j a (s) = (1 + s)~ 3 / 2 . In order to compare the above 
expressions to the contribution from the tail part of the 
probability densities Eqs. PfJ]) and (|Tfl) . we summarize 
the dependencies on t = t — to and to relative to the 
contribution from x > 1, Xq > 1 in Tab. (JTTJ) . The contri- 
butions for x < 1 or xq < 1 are negligible for long times 
r and to and Eq. indeed gives the leading order of 
the correlation function. 



Appendix B: Discussion of the functions f a (s) and 

In Eq. (jl5|) we have introduced the two functions f a (s) 
and g Q (s), 

fa{s) = J d yy 2 e~ r 1F1 Q,a + l,y 2 ^ T (a,y 2 s) , 

g a {s) = J™ dyy 2 e-y\¥ 1 (~, a + 1, y 2 ) e^ 8 . (Bl) 

We now want to find simpler expressions for these func- 
tions that are valid in the limit of small, respectively 
large, s and thus give us the limiting behavior of the cor- 
relation function Eqs. and (|T7)) . First, we divide the 
integral at y = 1 into, 

/a(s) - 3l(s)+ J 2 (S) 

(We-^Fi Q,a + l,y 2 )r(a l2/ 2 S ). 

(B2) 

For small s, the argument of the incomplete Gamma 
function in the first integral is small and we may ap- 
proximate T(a,y 2 s) ~ T(a), so that, 



J\yy 2 e-y 2 1 F 1 Q, a + l,/)r 



(a) , (B3) 



which is a constant independent of s. In the second inte- 
gral, we introduce the variable z = ^J~s y to get, 

3 f°° 2 ^ /3 z 2 \ 

3 2 {s) = s~2 J ^ dzz 2 e~^ 1F1 \^-,a + 1, —J r(a,z 2 ). 



Since the argument of the hypergeometric function is 
large except for very small z (which gives us another 



constant contribution), we may use the large- argument 
expansion [Hj], Eq. (13.5.1) 

lPl (!, a+iy ) K S-+V'' V - (B5) 



/7T 



to write 



J 2 (s) 



2r(a + l) __ s 



dzz d -^r(a,z^) . (B6) 



This integral can be calculated and yields in the limit of 
small s, 



Ms) 



r(a+ 1) a _ 2 , 2 -a\ 

-s ( 1 — 1 [pt)s J 



>/5F(2-a)' 



(B7) 



The limiting behavior of this expression now depends on 
the value of a: For a > 2 the second term dominates for 
small s and ^2(5) is constant, while for a < 2, the first 
term dominates and we have ^2(5) oc s a ~ 2 . For small s, 
we thus have, 

/a(^{ rr (B8) 

r(a+l) nCC -2 



for a < 2 . 



The value of the constant for a > 2 is obtained by eval- 
uating / a (0) using Mathematica. 

For large s, we introduce the variable z — \fsy in 
3i(s), 

_i /"^* * 2 ^ /3 z 2 \ 

3i(s)=s~ 2 / dzz 2 e » 1F1 ( -,a + l, — Jr(o;,z 2 ). 

(B9) 

The argument of the exponential and hypergeometric 
functions now is small except for very large z, for which 
the incomplete Gamma function is exponentially small 
(see below). For small arguments, both the exponential 
and hypergeometric functions are approximately 1, so we 
have, 

rs/i 

3i(s)~s-' dzz 2 T(a,z 2 ), (BIO) 
J 



(B4) which for large s reduces to 



3i( s ) o s 



(Bll) 
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In ^2(3), we expand the incomplete Gamma function for 
large arguments ([35j, Eq. (6.5.32)) 



T(a,y)~y a - 1 e-y 
and use Eq. (|B5|) to get, 

Ms) - ^±^.s- 
Evaluating the integral gives, 

3 2 ( S ) * £^±ll s «- 2e - , 



(B12) 



(B13) 



(B14) 



which vanishes exponentially for large s. So, for large s, 
we have, 



Us) 



r a 



-s 2 



(B15) 



From Eqs. (|B8l) and (IB15|) . we have for /„(s) 



^ r( 4r^-1) ( "" 2) for s « 1 and a > 2 
0=^ ^^)^" 22 for S «landa<2 
r(g+f) 



-s 2 



In a similar manner, we obtain, 



7o(«) 



r(q+l)r(2-q) „_i 



_ 1 

S 2" 



for s ^ 1 



for s <C 1 
for s > 1 



(B16) 



(B17) 



Appendix C: Discussion of the integral over the 
correlation function in Eq. (|54[) 

Before evaluating the double time integral over the cor- 
relation function C(t", t') in Eq. ([5"4")) . we have to address 
the problem that there are regions where r = t" — t! is 
not large and thus the approximate expression Eq. (|4"4")l 
is not valid. In these regions, we make use of the fact 
that C(t",t') < {x{t') 2 ) for t" > t' to provide an upper 
bound on the value of the correlation function. We write, 



dt" f dt'C(t",t') 
t* Jt* 

t pt"— A 

dt" / dt'C(t",t') 

t*+A Jt* 



= 2 



t ft" 

dt" / dt'C(t",t') 

f Jt"-A 
t*+A ft* 

dt" / dt'C(t",t') 

t* Jt"-A 



(CI) 



The hrst of the three integrals on the right-hand side now 
satisfies r = t" — t' > A with A chosen such that the ap- 
proximation (I44p of the correlation function is accurate. 



The time scale A can in principle be obtained by com- 
paring the result (|44p with the exact correlation function 
(obtained e.g. numerically). However, the asymptotic 
behavior will turn out to be independent of A in most 
cases. In the remaining two integrals, we use the esti- 
mate C(t",t') < (x(t') 2 ) and the expression fl32]) for the 
second moment to obtain, 



t ft" ft*+A ft* 

dt" / dt' - dt" / dt' 

t* Jt"—A Jt* Jt"-A 



C(t",t') 



3a— 6 
, , 1 



Zr(a)(2-a)(3-a) 
2 



i=^(4Dt)tA 



for a > 2 
— (4Dt)- a+2 tA forl<a<2 
for a < 1 



(C2) 



in the limit of large t (specifically t ^> t* and t ^> A). 
Using the general form of the correlation function (|44|) . 



t" - 1' 
t> 



c(t",t') ~ (t" -t'y<j> 



we can express the first integral in Eq. (IC1|) as 



t ft"-A 
dt" / dt'C{t",t') 
Jt* 



(C3) 



(C4) 



Making the change of variable s = (t" — t')/t' in the 
i'-integral, we further have 



dt"t"^ 



(s + l) 



M+2 



000 



(C5) 



where the function (j>(s) is given by Eq. (|45l) . 
/ oo «iW 2 e-^ 1 Fi(| l a + l, l r 1 ) 



xr(a,y^s) 
J^dyy 2 1 F 1 (la + l,y 2 ) 



for a > 1 



xe 



V(«+i) 



for a < 1 



(C6) 



with the asymptotic behavior (see appendix IB1 . 

for s 3> 1 and a > 1 



0(a) 



s 2 

g-^-a f or s ^> 1 an d a < l 

const, for s <C 1 and a > 2 

s q-2 for s < 1 and 1 < a < 2 

,s _1 for s«l and a < 1 . 



(C7) 



Note that for now we omit any s-independent prefactors 
of 4>(s), these have been included in the result in section 
fVl The exponent [i is given by, 



2 — a for a > 1 
1 for a < 1 



(C8) 



14 



When t" is close to the lower boundary, the s-integral 
vanishes since the lower and upper boundary are the 
same. For this reason the main contribution comes from 
t" close to the upper boundary, that is t" s» t. Using 
t 3> A and t S> t* , we have for the s-integral 



ds ■ 



f 4 (a + l)"+ 2 

We split this integral at s = 1 into, 



(C9) 



In the limit t ^> to, we then get for the expression 
Eq. EU, 



dt"t"» +1 



Jto+A 






~ < 


i 4 




^ 3 



*0 



t" - A 

for a > 3 
i 4 ~" for 1< a < 3 
for a < 1 . 



(s + 1)^+ 2 



(C15) 



(is ■ 



S+ l)^ 2 



tf>(s) 



ds 



+ 



M+2 



In order to obtain the prefactor to these asymptotic 
Ms) . forms, we take the limit t" — > oo in the s-integral in 

Eq. (|C15|) for a < 3, since the integral is constant in this 
(CIO) limit: 



In the first integral, we use the small argument expansion 
of <j)(s), Eq. (|C7| and obtain, 



ds - — — < 

a (s + \y+ 2 



l-(A) 3 -° fora>2 



t" 



(Gil) 



for a < 2 



For large t", the leading contribution for the case a > 2 
depends on a: 



ds ■ 



T<K*) 



a. ( S + iy+ 2 

(A) 3 "" fora>3 

const. for 2 < a < 3 

const. for a < 2 . 



(C12) 



In the second integral in Eq. (IC10I) . we use the large 
argument expansion of <fi(s), 



as ■ 



1 (S + 1)^+2 

l-(£) _f fora>l 
l-(g)~ i_Q fora<l 

So, for large i", we have 



(C13) 



(s + 1) 



M+2 



const. 



(C14) 



1 f°° s^ 
M + Jo dS (s + l)"+ 2 



0(s) (C16) 



For a > 3, the s-integral grows as t" 3 as the lower bound- 
ary approaches zero (see Eq. (IC12p ) and we may not take 
the lower boundary to zero. However, since the main 
contribution stems from small values of s, we may ap- 
proximate <f>(s) rs (j)(0) and have, 



ds 



(s + l) 



M+2 



0(0) 1 



A / dt"t'^ +1 I V 



-At - 1 t 2 
-M-l 

Summarizing these results, we have 



(C17) 



r £M (4£)A) 2 - Q ^ 



V A/ 
for a > 3 



zr 



(a+l)r(a)(4-a) ( 4jD< ) 2 "Jo rfs (s+1)*-" /"( S ) 

for 1 < a < 3 



3r(a)rfi- Q ) 4 -° f Jo ds O+i) 3 9a ( s ) 
for a < 1 , 



(C18) 



where we use the subscript I to denote the contribution 
from the long-time behavior of the correlation function. 
The functions f a (s) and g a (s) are the ones defined in 
Eq. (|45l) . The other contribution comes from the short- 
time correlation function, which we have approximated 
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by a constant in Eq. (|C2[) . 

( 2a-l f±^~ 1 



3q-6 



for a > 2 



zr(a)(2-a)(3-a) 

x(4£>A)(4D^) 1 - Q forl<a<2 



l-a 
2 



(4DA) 



for a < 1 . 



(C19) 



Comparing the long and short time contributions, 
Eqs. (|C18[) and (|C19p . we see that the long-time con- 
tribution dominates for a < 3, while both contributions 
are of the same order for a > 3. As a consequence, we 
have, 



! (*)) 



for a > 3 



zr 



(a+l)r(a)(4-Q) ( 4 - D *) 2 "/o ^ S (s+l) 4 -° /"( S ) 

for 1 < a < 3 



3r(a)r(i-a) 4-0* Jo ds ( s+1 yt g a (s) 
for a < 1 . 



(C20) 



The prefactor c a for a > 3 cannot be obtained within 
this asymptotic analysis, since it will in general depend 
on the shape of the potential U{x) in the whole space in- 
stead of just the logarithmic U (x) ~ Uo In (a;) asymptotic 
large x behavior. c a can, however, be obtained from the 
Boltzmann equilibrium distribution Q. 
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